RNA Secondary Structure Prediction Research Articles

The potential up-regulation risk of 3′ UTR SNP (rs10787760 G > A) for the VAX1 gene is associated with NSCLP in the northwest Chinese population

AimsTo investigate the association between single nucleotide polymorphisms (SNPs) in 3′UTR region of VAX1, SYT14 and PAX7 genes and the risk of non-syndromic cleft palate (NSCLP) in a northwest Chinese population. Main methodsA case-control study was conducted in 406 normal controls and 399 NSCLP patients. Using iMLDRTM genotyping technology, eight SNPs of three genes ((rs10787760, rs7086344 at VAX1), (rs1010113, rs851114, and rs485874 at PAX7), and (rs61820397, rs4609425, rs12133399 at SYT14)) were genotyped to investigate the differences in alleles and genotype distribution frequencies between NSCLP patients and healthy controls. RNA Folding Form software was used to predict RNA secondary structure and expression vectors were constructed to explore the function of the relevant SNP. The effect of SNP polymorphism of gene transcription and translation was assessed using qPCR and Western blot analysis.Key findings.Among the eight SNPs of three genes, rs10787760 of VAX1 gene was found to be associated with an increased risk of NSCLP (OR = 1.341,CI = 1.004–1.790) and the GA genotype of rs10787760 increased the risk of cleft lip and/or palate (CL/P) about 1.42 times (p < 0.05), and carrying the A allele might increase the risk of NSCL/P in male (OR = 1.356, 95 % CI = 1.010–1.823). But there was no association observed with cleft palate only (CPO). Cell function experiments revealed that the G to A mutation in rs10787760 up-regulated GFP-VAX1 transcriptional level by 2.39 and 3.13 times in two cell lines respectively, and enhance the protein expression of the VAX1 gene further. RNA secondary structure study showed that the rs10787760 (G > A) had two different secondary structures in 3′UTR region. SignificanceThe rs10787760 variant in the 3′UTR region of VAX1 gene is associated with CL/P in northwest Chinese population. We hypothesize that the machanism of it might be caused by the RNA differenct fold in the 3′UTR region caused by the polymorphism of the gene. Level of evidenceOriginal Reports.

Accurate prediction of RNA secondary structure including pseudoknots through solving minimum-cost flow with learned potentials

Pseudoknots are key structure motifs of RNA and pseudoknotted RNAs play important roles in a variety of biological processes. Here, we present KnotFold, an accurate approach to the prediction of RNA secondary structure including pseudoknots. The key elements of KnotFold include a learned potential function and a minimum-cost flow algorithm to find the secondary structure with the lowest potential. KnotFold learns the potential from the RNAs with known structures using an attention-based neural network, thus avoiding the inaccuracy of hand-crafted energy functions. The specially designed minimum-cost flow algorithm used by KnotFold considers all possible combinations of base pairs and selects from them the optimal combination. The algorithm breaks the restriction of nested base pairs required by the widely used dynamic programming algorithms, thus enabling the identification of pseudoknots. Using 1,009 pseudoknotted RNAs as representatives, we demonstrate the successful application of KnotFold in predicting RNA secondary structures including pseudoknots with accuracy higher than the state-of-the-art approaches. We anticipate that KnotFold, with its superior accuracy, will greatly facilitate the understanding of RNA structures and functionalities.

Structural Impact of the Interaction of the Influenza A Virus Nucleoprotein with Genomic RNA Segments.

Influenza A viruses (IAVs) possess a segmented genome consisting of eight viral RNAs (vRNAs) associated with multiple copies of viral nucleoprotein (NP) and a viral polymerase complex. Despite the crucial role of RNA structure in IAV replication, the impact of NP binding on vRNA structure is not well understood. In this study, we employed SHAPE chemical probing to compare the structure of NS and M vRNAs of WSN IAV in various states: before the addition of NP, in complex with NP, and after the removal of NP. Comparison of the RNA structures before the addition of NP and after its removal reveals that NP, while introducing limited changes, remodels local structures in both vRNAs and long-range interactions in the NS vRNA, suggesting a potentially biologically relevant RNA chaperone activity. In contrast, NP significantly alters the structure of vRNAs in vRNA/NP complexes, though incorporating experimental data into RNA secondary structure prediction proved challenging. Finally, our results suggest that NP not only binds single-stranded RNA but also helices with interruptions, such as bulges or small internal loops, with a preference for G-poor and C/U-rich regions.

SparseRNAfolD: optimized sparse RNA pseudoknot-free folding with dangle consideration

MotivationComputational RNA secondary structure prediction by free energy minimization is indispensable for analyzing structural RNAs and their interactions. These methods find the structure with the minimum free energy (MFE) among exponentially many possible structures and have a restrictive time and space complexity (O(n3)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$O(n^3)$$\\end{document} time and O(n2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$O(n^2)$$\\end{document} space for pseudoknot-free structures) for longer RNA sequences. Furthermore, accurate free energy calculations, including dangle contributions can be difficult and costly to implement, particularly when optimizing for time and space requirements.ResultsHere we introduce a fast and efficient sparsified MFE pseudoknot-free structure prediction algorithm, SparseRNAFolD, that utilizes an accurate energy model that accounts for dangle contributions. While the sparsification technique was previously employed to improve the time and space complexity of a pseudoknot-free structure prediction method with a realistic energy model, SparseMFEFold, it was not extended to include dangle contributions due to the complexity of computation. This may come at the cost of prediction accuracy. In this work, we compare three different sparsified implementations for dangle contributions and provide pros and cons of each method. As well, we compare our algorithm to LinearFold, a linear time and space algorithm, where we find that in practice, SparseRNAFolD has lower memory consumption across all lengths of sequence and a faster time for lengths up to 1000 bases.ConclusionOur SparseRNAFolD algorithm is an MFE-based algorithm that guarantees optimality of result and employs the most general energy model, including dangle contributions. We provide a basis for applying dangles to sparsified recursion in a pseudoknot-free model that has the potential to be extended to pseudoknots.

Prediction of RNA Secondary Structure Using Butterfly Optimization Algorithm

AbstractRibonucleic acid (RNA) structure is vital to its ability to function within the cell. The ability to predict RNA structure is essential to implementing new medications and understanding genetic illnesses. It is also important in synthetic and computational biology. All these functions are directly related to its secondary structure. Also prediction of RNA secondary structure process is the most significant step to determining the tertiary structure of RNA. On account of this, prediction of secondary structure of RNA is the crying topic in bioinformatics. In this research, we present the swarm-based metaheuristic Butterfly Optimization Algorithm (BOA) method for predicting the secondary structure of RNA. The main feather of the BOA is that it can conduct both local and global search simultaneously. According to the problem perspective, we have redesigned the operators of BOA to perform global and local search operations in different ways. We have followed a thermodynamic model for the selection of the stable secondary structure with minimum Gibbs free energy. Predicting the minimum free energy value we also developed an “Optimize” function to search the new optimize structure. This function increases the prediction efficiency, creating new stable structure and also decreases the time complexity of global searching procedure. We have used a public dataset to perform the prediction operation. To accuse our prediction efficiency, we have compared our outcomes to existing popular algorithms. The result shows that the proposed approach can predict secondary RNA structure better than other state-of-the-art algorithms.

NNDB: An Expanded Database of Nearest Neighbor Parameters for Predicting Stability of Nucleic Acid Secondary Structures

Nearest neighbor thermodynamic parameters are widely used for RNA and DNA secondary structure prediction and to model thermodynamic ensembles of secondary structures. The Nearest Neighbor Database (NNDB) is a freely available web resource (https://rna.urmc.rochester.edu/NNDB) that provides the functional forms, parameter values, and example calculations. The NNDB provides the 1999 and 2004 set of RNA folding nearest neighbor parameters. We expanded the database to include a set of DNA parameters and a set of RNA parameters that includes m6A in addition to the canonical RNA nucleobases. The site was redesigned using the Quarto open-source publishing system. A downloadable PDF version of the complete resource and downloadable sets of nearest neighbor parameters are available.

Prediction of the effects of the top 10 synonymous mutations from 26645 SARS-CoV-2 genomes

Background The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) had led to a global pandemic since December 2019. SARS-CoV-2 is a single-stranded RNA virus, which mutates at a higher rate. Multiple works had been done to study nonsynonymous mutations, which change protein sequences. However, there is little study on the effects of SARS-CoV-2 synonymous mutations, which may affect viral fitness. This study aims to predict the effect of synonymous mutations on the SARS-CoV-2 genome. Methods A total of 26645 SARS-CoV-2 genomic sequences retrieved from Global Initiative on Sharing all Influenza Data (GISAID) database were aligned using MAFFT. Then, the mutations and their respective frequency were identified. Multiple RNA secondary structures prediction tools, namely RNAfold, IPknot++ and MXfold2 were applied to predict the effect of the mutations on RNA secondary structure and their base pair probabilities was estimated using MutaRNA. Relative synonymous codon usage (RSCU) analysis was also performed to measure the codon usage bias (CUB) of SARS-CoV-2. Results A total of 150 synonymous mutations were identified. The synonymous mutation identified with the highest frequency is C3037U mutation in the nsp3 of ORF1a. Of these top 10 highest frequency synonymous mutations, C913U, C3037U, U16176C and C18877U mutants show pronounced changes between wild type and mutant in all 3 RNA secondary structure prediction tools, suggesting these mutations may have some biological impact on viral fitness. These four mutations show changes in base pair probabilities. All mutations except U16176C change the codon to a more preferred codon, which may result in higher translation efficiency. Conclusion Synonymous mutations in SARS-CoV-2 genome may affect RNA secondary structure, changing base pair probabilities and possibly resulting in a higher translation rate. However, lab experiments are required to validate the results obtained from prediction analysis.

Limits of experimental evidence in RNA secondary structure prediction.

Limits of experimental evidence in RNA secondary structure prediction.

Concurrent prediction of RNA secondary structures with pseudoknots and local 3D motifs in an integer programming framework.

The prediction of RNA structure canonical base pairs from a single sequence, especially pseudoknotted ones, remains challenging in a thermodynamic models that approximates the energy of the local 3D motifs joining canonical stems. It has become more and more apparent in recent years that the structural motifs in the loops, composed of noncanonical interactions, are essential for the final shape of the molecule enabling its multiple functions. Our capacity to predict accurate 3D structures is also limited when it comes to the organization of the large intricate network of interactions that form inside those loops. We previously developed the integer programming framework RNA Motifs over Integer Programming (RNAMoIP) to reconcile RNA secondary structure and local 3D motif information available in databases. We further develop our model to now simultaneously predict the canonical base pairs (with pseudoknots) from base pair probability matrices with or without alignment. We benchmarked our new method over the all nonredundant RNAs below 150 nucleotides. We show that the joined prediction of canonical base pairs structure and local conserved motifs (i) improves the ratio of well-predicted interactions in the secondary structure, (ii) predicts well canonical and Wobble pairs at the location where motifs are inserted, (iii) is greatly improved with evolutionary information, and (iv) noncanonical motifs at kink-turn locations. The source code of the framework is available at https://gitlab.info.uqam.ca/cbe/RNAMoIP and an interactive web server at https://rnamoip.cbe.uqam.ca/.

New approximation algorithms for RNA secondary structures prediction problems by local search

This paper investigates two combinatorial problems from RNA secondary structure prediction with arbitrary pseudoknots. Given a RNA sequence and a set of base pairs, two parallel and adjacent base pairs constitute a stacking. The Maximum Stacking Base Pairs problem (MSBP) aims at finding a maximum number of based pairs, all of which form stackings, while the Maximum Base Pair Stackings problem (MBPS) is to find a maximum number of stackings. Both problems are NP-hard. We present two new approximation algorithms for the two problems by local search methods. For MSBP, the approximation factor is improved from 52 to 73; as for the MBPS, the approximation factor is improved from 83 to 52.

Competitive Influence of Alkali Metals in the Ion Atmosphere on Nucleic Acid Duplex Stability.

The nonspecific atmosphere around nucleic acids, often termed the ion atmosphere, encompasses a collection of weak ion-nucleic acid interactions. Although nonspecific, the ion atmosphere has been shown to influence nucleic acid folding and structural stability. Studies investigating the composition of the ion atmosphere have shown competitive occupancy of the atmosphere between metal ions in the same solution. Many studies have investigated single ion effects on nucleic acid secondary structure stability; however, no comprehensive studies have investigated how the competitive occupancy of mixed ions in the ion atmosphere influences nucleic acid secondary structure stability. Here, six oligonucleotides were optically melted in buffers containing molar quantities, or mixtures, of either XCl (X = Li, K, Rb, or Cs) or NaCl. A correction factor was developed to better predict RNA duplex stability in solutions containing mixed XCl/NaCl. For solutions containing a 1:1 mixture of XCl/NaCl, one alkali metal chloride contributed more to duplex stability than the other. Overall, there was a 54% improvement in predictive capabilities with the correction factor compared with the standard 1.0 M NaCl nearest-neighbor models. This correction factor can be used in models to better predict RNA secondary structure in solutions containing mixed XCl/NaCl.

Accelerating prediction of RNA secondary structure using parallelization on multicore architecture

Accelerating prediction of RNA secondary structure using parallelization on multicore architecture

Automated design of dynamic programming schemes for RNA folding with pseudoknots

Although RNA secondary structure prediction is a textbook application of dynamic programming (DP) and routine task in RNA structure analysis, it remains challenging whenever pseudoknots come into play. Since the prediction of pseudoknotted structures by minimizing (realistically modelled) energy is NP-hard, specialized algorithms have been proposed for restricted conformation classes that capture the most frequently observed configurations. To achieve good performance, these methods rely on specific and carefully hand-crafted DP schemes. In contrast, we generalize and fully automatize the design of DP pseudoknot prediction algorithms. For this purpose, we formalize the problem of designing DP algorithms for an (infinite) class of conformations, modeled by (a finite number of) fatgraphs, and automatically build DP schemes minimizing their algorithmic complexity. We propose an algorithm for the problem, based on the tree-decomposition of a well-chosen representative structure, which we simplify and reinterpret as a DP scheme. The algorithm is fixed-parameter tractable for the treewidth tw of the fatgraph, and its output represents a O(ntw+1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathcal {O}}(n^{tw+1})$$\\end{document} algorithm (and even possibly O(ntw)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathcal {O}}(n^{tw})$$\\end{document} in simple energy models) for predicting the MFE folding of an RNA of length n. We demonstrate, for the most common pseudoknot classes, that our automatically generated algorithms achieve the same complexities as reported in the literature for hand-crafted schemes. Our framework supports general energy models, partition function computations, recursive substructures and partial folding, and could pave the way for algebraic dynamic programming beyond the context-free case.

RNAknot2S: a New RNA Secondary Structure Prediction Method that Combines the Comparative Approach with Genetic Algorithms

RNAknot2S: a New RNA Secondary Structure Prediction Method that Combines the Comparative Approach with Genetic Algorithms

AliNA - a deep learning program for RNA secondary structure prediction.

Nowadays there are numerous discovered natural RNA variations participating in different cellular processes and artificial RNA, e.g., aptamers, riboswitches. One of the required tasks in the investigation of their functions and mechanism of influence on cells and interaction with targets is the prediction of RNA secondary structures. The classic thermodynamic-based prediction algorithms do not consider the specificity of biological folding and deep learning methods that were designed to resolve this issue suffer from homology-based methods problems. Herein, we present a method for RNA secondary structure prediction based on deep learning - AliNA (ALIgned Nucleic Acids). Our method successfully predicts secondary structures for non-homologous to train-data RNA families thanks to usage of the data augmentation techniques. Augmentation extends existing datasets with easily-accessible simulated data. The proposed method shows a high quality of prediction across different benchmarks including pseudoknots. The method is available on GitHub for free (https://github.com/Arty40m/AliNA).

In vivo-like nearest neighbor parameters improve prediction of fractional RNA base-pairing in cells.

We conducted a thermodynamic analysis of RNA stability in Eco80 artificial cytoplasm, which mimics in vivo conditions, and compared it to transcriptome-wide probing of mRNA. Eco80 contains 80% of Escherichia coli metabolites, with biological concentrations of metal ions, including 2 mM free Mg2+ and 29 mM metabolite-chelated Mg2+. Fluorescence-detected binding isotherms (FDBI) were used to conduct a thermodynamic analysis of 24 RNA helices and found that these helices, which have an average stability of -12.3 kcal/mol, are less stable by ΔΔGo37 ∼1 kcal/mol. The FDBI data was used to determine a set of Watson-Crick free energy nearestneighbor parameters (NNPs), which revealed that Eco80 reduces the stability of three NNPs. This information was used to adjust the NN model using the RNAstructure package. The in vivo-like adjustments have minimal effects on the prediction of RNA secondary structures determined in vitro and in silico, but markedly improve prediction of fractional RNA base pairing in E. coli, as benchmarked with our in vivo DMS and EDC RNA chemical probing data. In summary, our thermodynamic and chemical probing analyses of RNA helices indicate that RNA secondary structures are less stable in cells than in artificially stable in vitro buffer conditions.

Automatic recognition of complementary strands: lessons regarding machine learning abilities in RNA folding.

Introduction: Prediction of RNA secondary structure from single sequences still needs substantial improvements. The application of machine learning (ML) to this problem has become increasingly popular. However, ML algorithms are prone to overfitting, limiting the ability to learn more about the inherent mechanisms governing RNA folding. It is natural to use high-capacity models when solving such a difficult task, but poor generalization is expected when too few examples are available. Methods: Here, we report the relation between capacity and performance on a fundamental related problem: determining whether two sequences are fully complementary. Our analysis focused on the impact of model architecture and capacity as well as dataset size and nature on classification accuracy. Results: We observed that low-capacity models are better suited for learning with mislabelled training examples, while large capacities improve the ability to generalize to structurally dissimilar data. It turns out that neural networks struggle to grasp the fundamental concept of base complementarity, especially in lengthwise extrapolation context. Discussion: Given a more complex task like RNA folding, it comes as no surprise that the scarcity of useable examples hurdles the applicability of machine learning techniques to this field.

New evidence of consistency between phylogeny and morphology for two taxa in ciliated protists, the subclasses Oligotrichia and Choreotrichia (Protista, Ciliophora)

Marine planktonic ciliates are largely oligotrichs and choreotrichs, which are two subclasses of the class Spirotrichea. The current phylogenetic assignments of oligotrichs and choreotrichs are inconsistent with previous results based on morphological features, probably hindered by the limited information from a single gene locus. Here we provide 53 new sequences from small subunit ribosomal RNA (SSU rDNA), ITS1-5.8S rDNA-ITS2, and large subunit ribosomal RNA (LSU rDNA) gene loci in 25 oligotrich and choreotrich species. We also predict RNA secondary structures for the ITS2 regions in 55 species, 48 species of which are reported for the first time. Based on these novel data, we make a more comprehensive phylogenetic reconstruction, revealing consistency between morphological taxonomy and an updated phylogenetic system for oligotrichs and choreotrichs. With the addition of data from ciliature patterns and genes, the phylogenetic analysis of the subclass Oligotrichia suggests three evolutionary trajectories, among which: 1) Novistrombidium asserts an ancestral ciliary pattern in Oligotrichia; 2) the subgenera division of Novistrombidium and Parallelostrombidium are fully supported; 3) the three families (Tontoniidae, Pelagostrombidiidae and Cyrtostrombidiidae) all evolved from the most diverse family Strombidiidae, which explains why strombidiids consistently form polyphyletic clades. In the subclass Choreotrichia, Strombidinopsis likely possesses an ancestral position to other choreotrichs, and both phylogenetic analysis and RNA secondary structure prediction support the hypothesis that tintinnids may have evolved from Strombidinopsis. The results presented here offer an updated hypothesis for the evolutionary history of oligotrichs and choreotrichs based on new evidence obtained by expanding sampling of molecular information across multiple gene loci.

Deciphering phenotyping, DNA barcoding, and RNA secondary structure predictions in eggplant wild relatives provide insights for their future breeding strategies

Eggplant or aubergine (Solanum melongena L.) and its wild cousins, comprising 13 clades with 1500 species, have an unprecedented demand across the globe. Cultivated eggplant has a narrow molecular diversity that hinders eggplant breeding advancements. Wild eggplants need resurgent attention to broaden eggplant breeding resources. In this study, we emphasized phenotypic and genotypic discriminations among 13 eggplant species deploying chloroplast–plastid (Kim matK) and nuclear (ITS2) short gene sequences (400–800 bp) at DNA barcode region followed by ITS2 secondary structure predictions. The identification efficiency at the Kim matK region was higher (99–100%) than in the ITS2 region (80–90%). The eggplant species showed 13 unique secondary structures with a central ring with various helical orientations. Principal component analysis (PCoA) provides the descriptor–wise phenotypic clustering, which is essential for trait–specific breeding. Groups I and IV are categorized under scarlet complexes S. aethiopicum, S. trilobatum, and S. melongena (wild and cultivated). Group II represented the gboma clade (S. macrocarpon, S. wrightii, S. sisymbriifolium, and S. aculeatissimum), and group III includes S. mammosum, and S. torvum with unique fruit shape and size. The present study would be helpful in genetic discrimination, biodiversity conservation, and the safe utilization of wild eggplants.

Mono-valent salt corrections for RNA secondary structures in the ViennaRNA package

BackgroundRNA features a highly negatively charged phosphate backbone that attracts a cloud of counter-ions that reduce the electrostatic repulsion in a concentration dependent manner. Ion concentrations thus have a large influence on folding and stability of RNA structures. Despite their well-documented effects, salt effects are not handled consistently by currently available secondary structure prediction algorithms. Combining Debye-Hückel potentials for line charges and Manning’s counter-ion condensation theory, Einert et al. (Biophys J 100: 2745-2753, 2011) modeled the energetic contributions of monovalent cations on loops and helices.ResultsThe model of Einert et al. is adapted to match the structure of the dynamic programming recursion of RNA secondary structure prediction algorithms. An empirical term describing the salt dependence of the duplex initiation energy is added to improve co-folding predictions for two or more RNA strands. The slightly modified model is implemented in the ViennaRNA package in such way that only the energy parameters but not the algorithmic structure is affected. A comparison with data from the literature show that predicted free energies and melting temperatures are in reasonable agreement with experiments.ConclusionThe new feature in the ViennaRNA package makes it possible to study effects of salt concentrations on RNA folding in a systematic manner. Strictly speaking, the model pertains only to mono-valent cations, and thus covers the most important parameter, i.e., the NaCl concentration. It remains a question for future research to what extent unspecific effects of bi- and tri-valent cations can be approximated in a similar manner.AvailabilityCorrections for the concentration of monovalent cations are available in the ViennaRNA package starting from version 2.6.0.